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PACS. 68.35.Ja - Surface and interface dynamics and vibrations. 
I . PACS. 81.10.Blc- Growth from vapor. 

' ^ ' PACS. 64.60.Ht - Dynamic critical phenomena. 

o 

»^ ' Abstract. - We study, both with numerical simulations and theoretical methods, a cellular 

►^ , automata model for continuum equations describing growth processes in the presence of an 

external flux of particles. As a result of local instabilities we find a fractal regime of growth for 
small external fluxes. The growing tip is selected with probability proportional to the curvature 
ivj ' in the point. A parameter p gives the probability of lateral growth on the tip. The value of p 

^y. ^ determines the fractal dimension of the aggregate. Furthermore, for each value of p a cross-over 

J.— s , between two different fractal dimensions is observed. Instead, the roughness exponent x oi the 

^— «, ' aggregates does not depend on p (x — 0.5). Fixed scale transformation approach is applied to 

(<-^ ' compute theoretically the fractal dimension for one of the branches of the structure. 

B 

M , From a technological point of view the growth by vapor atom deposition is the most common 

O ' way to produce films jl| . It is known that local instabilities play a decisive role in the formation 

of patterns observed in experiments, such as grooves, mounds and metallurgical junctions H. 
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These kinds of phenomena are usually described by the following continuum equation 



^- hix)^J~^V'hix), (1) 

C^ , 

where h{x) represents the surface height corresponding to the substrate coordinate x; J is 
the incident flux of particles. The Laplacian term with the negative sign plays the role of an 
antidilFusion, and it is responsible for the instability. The precise meaning of £, depends on the 
particular kind of instability involved in the process. The J term in eq.(E|) can be ruled out 
with a change of coordinates in a frame comoving with the surface. Then h can be considered 
just as the fluctuation around the average. The resulting equation is 

h{x,t) = -V^h{x,t). (2) 
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The solution of (Q) in Fourier space is 

h{k,t)^h{k,0)e^^'*. (3) 

Therefore all the modes present at the origin are exponentially amplified. In order to control 
this instability, people usually introduce a term proportional to \7'^h{x) in the eq. (m). This 
brings to 

/i(fc, t) = iJ?h{k, t) - \k^h{k, t) (4) 

whose solution now is 

h(k,t)^h{k,Q)e^''^-^^\. (5) 

Therefore only the long wavelength modes are still unstable. The result of this process is the 
formation of columns separated by very deep grooves. 

An alternative continuum approach is represented by the Reparametrization Invariant 
Formulation of the surface growth equation Q . The basic idea is that the growth of surfaces 
only depends on its intrinsic properties. Therefore eq.(|l|) can be rewritten as 

'r = h{J^^H) (6) 

where r is the position vector of the surface, h the local normal unit vector and H the local 
curvature. J and f are defined as before. Within this scheme the 'V'/i term is substituted by 
the Beltrami-Laplace operator for the curvature. 

When the tips of two columns join together the deep groove below them disappears and 
the effective resulting surface is almost flat|5|]. In this spirit, a simpler way to approach the 
study of surface growth is represented by the introduction of lattice models able to keep the 
relevant properties of the phenomenon. Many models have been proposed in the past years: 
the Eden model [g|, ballistic deposition models 0, Solid-On-Solid (SOS) models Q and others 
(for a review on surface growth models see also Ref.|^, 10 1). 



Here, we present a lattice model to verify the role of the surface curvature in these growth 
processes. Consider a surface subject to an incoming flux of particles. The model is based 
on the following two rules. Firstly, the probability that a particle stops in a given site is 
proportional to the local negative curvature of the surface. For each point the curvature is 
computed in a first approximation as the Laplacian function in the point. In this way, the 
growth probability at t on the tip af, t is given by 

P(x, t) oc [2h{x, t) - h(x - 1, t) - h[x + 1, t)] (7) 

Secondly, a particle sticks right on the tip with probability p while with probability 1 — p 
it can stick laterally. We introduce in this way the possibility to develop overhangs. Under 
this respect this model can be considered a generalization of the growth model proposed in 
Ref.||ll[ for sticky particles. The surface height h{i) on a site i corresponds to the position of 
the highest particle on that column. We will consider here the limit of vanishing external flux, 
which is mimicked allowing just one particle to stick at a time. 

This model of aggregation does not give rise to compact structures, but instead to fractal 
ones merging themselves to form a continuous surface. In the system, every now and then, a tip 
whose value of the Laplacian operator becomes larger and larger, spans uniformly the system 
size. When the growth reaches a boundary there is a new flat surface where the deposition 
proceeds. 

To reduce finite size effects, in our simulations we considered periodic boundary conditions 
along the horizontal direction of the system. In Fig.y we present some typical patterns for 
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size L X 2L {L = 256) and for different values of p. There is a transient regime, that we cah 
the scaling regime in analogy with the dielectric breakdown model (DBM) |12| . Here several 
branches compete with each other in the growth process. After that, in the steady state, a single 
branch prevails and grows forever. Due to the screening induced by the shadowing between 
competing branches, the system develops fractal properties in the scaling regime where we 
focus our studies. We have performed several numerical simulations of this model for different 
values of p and for system sizes up to L = 512, and height 2L. The results presented here refer 
to average on about 100 realization for every value of p and L = 256, 512. 

In this model the sample grows as a compact set of different columns as observed experi- 
mentally in [O, nj]. Authors, in |14|, conjecture an exponential distribution for the tip height 
from their experimental results on columnar structured gold electrodeposits. 

D{h) oc exp{-bh) (8) 

where b is the inverse of a characteristic length. The same shape of D{h) was previously 



observed in 1 13 , through scanning tunneling microscopy analysis of vapor deposited gold 
films. We computed D(h) for our model with p = I, L = 256 and maximum tip height 
hmax = 256, 512. The results reproduce the same exponential behavior of D(h) observed in 
[ [l3| [l4[ . The characteristic height b^^ turns out to be actually proportional to the maximum 

height hmax {b^^ = C~^hmax ~ ii^). 

In order to describe the clusters obtained we also focused on the roughness exponent x{p) of 
the surface. The roughness exponent is defined as the ensemble averaged width of the interface 

W{1, t) = ((;^(^, t) - {h{t, t)))Y/' ~ l^'fit'/yi) (9) 

where z is the dynamical exponent, the angular brackets denote the average over all portions of 
the interface of length I for a given realization of total length L, and over different realizations, 
and / is a scaling function such that f{y) ^ y^ ioi y « 1 and f{y) = const for y >> 1. The 
exponent (3 = x/z describes the transient roughening. In the asymptotic time limit on gets: 

W{l,t^oo) = W{l)(xl^ (10) 

Interestingly enough, the value of this quantity does not depend on the value of p as one can 
see from Table |. A plot of this scaling is shown in Fig.y. 

A dependence on p affects, instead, the fractal dimension of the aggregate. In our simula- 
tions we studied the lower strip where the growth process can be considered fully screened by 
the shadowing effect. In Fig.^ we show the plot of N{1), number of occupied boxes of size I, 
versus I in the scaling regime for different values of p. In Tab.| we report the measured values 
for these exponents. The result is that the value of Df strongly depends on the parameter p. 
Moreover, for all values of p a cross-over between two different fractal dimensions is observed. 
In particular the region of small I are characterized by a value of Dl"™ changing from 1.20(3) 
to 1.08(3) as p passes from 0.1 to 0.9. Larger regions are characterized by a similar trend with 
the value of DV^ ' changing from 1.30(3) to 1.12(3) as p passes from 0.1 to 0.9 

As regards the trend of this variation one can conclude that in the limit p ^ 1 the growth is 
dominated by the highest tip that grows forever. In this condition the growth of the aggregate 
leaves large empty regions that strongly reduce the fractal dimension of the aggregate. The 
possibility of lateral growth instead, delays this effect in such a way to form an intermediate 
region where one can study non trivial fractal properties. 

In order to obtain an analytical computation to reproduce the behavior of Z?l°™, we used 
a theoretical scheme of a novel type that can deal with irreversible dynamic processes in 
a reasonably systematic way. The Fixed Scale Transformation (FST) |13] focuses on the 
dynamics at a given scale and allows to compute the nearest neighbors correlations at this 
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scale by lattice diagrams. The structure that we want to study becomes well defined in the 
following limits: 

asymptotic time limit {t —>■ cxd): only regions of the system where the growth is frozen 
are considered; 

large scale limit (r -^ oo) to avoid the presence of finite size effects 

In the FST scheme, these two limits are reached through the fixed point convergence of the 
lattice diagrams after a certain order and considering the system at a coarse grained scale 
m&. In order to perform the first limit one defines an elementary cell of two sites. The two 
configurations in Fig.H are the only ones subject to growth, with an associated probability to 
have a transition from one to the other in the growth direction. The fixed point of the process 
(that can be regarded as the limit t -^ cx)) is obtained looking at the eigenvector corresponding 
to the eigenvalue 1 (which is also the largest one) 

where Ci and C2 are the probabilities of occurrence of configurations of the type 1 and 2 
respectively. 

The average number of boxes occupied by the structure in a change of scale ^ — > ^ is 

< n >= Ci + 2C2. (12) 

Since we are studying the intersection set between the structure and a line defining the growth 
front, the fractal dimension of the structure is simply 

Df = l + Dt (13) 

where Dt is the transverse fractal dimension related to < n > by the formula: 

r, _ ln<n> _ ln(Ci+2C2) ,^., 

^* ~ n^7(2r ~ — M^^ — ■ ^ ^ 

Given this theoretical framework, it is very easy to develop the calculations for the present 
case. Consider the graph in FigJj; as a first approximation of the boundary conditions, one can 
consider the open ones, where the pattern has no other occupied sites around. The opposite 
condition is represented by the closed ones in which the entire pattern is replied. In both 
cases the growth probabilities are normalized to one inside the growth column. Usually, a 
satisfactory approximation consists in considering only these two limiting cases, where the two 
fixed value of the configuration probabilities are given by: 



Aqi + 2Mi[-lM^f-^{lM°f-M^i-2Mi)^~4AMii 

^1 = 2A 

C2 = 1-Cr (15) 

where A = Affj + -^^21 ~ 3/2 (M"!" + M!^^). By using the scheme shown in Fig.|| one obtains 
the values of the matrix elements in the FST scheme for open and closed boundary conditions. 

M^f^i^ ; Mil = 1- Mil (16) 

^"^' = Yh ' ^^^2T = 1-M2T (17) 
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Mt^^\-p ■ Mf, = l-Mf, (18) 

Mi=p ; M^i,^l-M^i (19) 

Two sources of screening are present in our model: one is due to the fact that only the tips can 
grow, and this is a local effect which is accounted for in the FST computation. The other, more 
important, is shadowing, which is highly non-local, depending on the whole structure of the 
aggregate. This effect cannot be easily accounted in the FST scheme ||l5|. This reflects on the 
results of the computation, especially for low values of p, where shadowing is more important. 
Instead for as p — > 1 there is a reasonable agreement between the FST computation and 
numerical simulations (last column of Table Q) . 

To conclude, we have presented a simple cellular automata model, which describes fractal 
surface growth in the limit of vanishingly small external flux, under the effect of a local 
instability. The model shows non-trivial, interesting fractal and roughness properties, which 
are studied both numerically and analytically. We acknowledge useful discussion with R. Ball, 
P. De Los Rios, A. Maritan, L. Pietronero, and the support of the European Network contract 
FMRXCT980183. 
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Table I. - Fractal dimension for the lower and upper scaling region and roughness exponent for 
different values of p and system size L = 256, 512. 



p 


D\°^(p) 


i?P(p) 


x{p) 


Di°^FST{p) 


0.1 


1.20 ±0.03 


1.30 ±0.03 


0.54 ±0.03 


1.91 


0.33 


1.18 ±0.03 


1.28 ±0.03 


0.55 ±0.04 


1.66 


0.5 


1.17 ±0.03 


1.24 ±0.02 


0.50 ±0.02 


1.46 


0.7 


1.14 ±0.02 


1.20 ±0.02 


0.48 ± 0.04 


1.24 


0.9 


1.08 ±0.02 


1.12 ±0.02 


- 


1.08 




Fig. 1. - (a)-(d): Typical patterns for different values of p (the system size and the values of p are 
shown in the legends). 
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Fig. 2. - Scaling behavior of the squared width W'^if) width of portions of size / of the aggregates 
surface, defined as the path of the highest occupied sites for each column, for different p. We get a 
roughness exponent x around 1/2 for all values of p. 
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Fig. 3. - (a)-(d): Scaling behavior of the number of occupied boxes A'^(/) of size I versus /, for different 
values of p, horizontal size L = 256, 512 and height h = 2L. The box counting is made on a strip of 
height h = L/2. 
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Fig. 4. - Lattice path integrals for the two types of boundary conditions considered. Please note that 
both in OBC than in the CBC, all the processes starting from configuration of the type C2 gives the 
same results of processes starting from configuration Ci. In fact, if one consider the site b black, it 
is impossible to have growth process starting from it, since the Laplacian operator in the point is 
negative 



